Measurement of the forbidden Gs^^Sq — > 5d6s ^Z?i magnetic-dipole transition 

amplitude in atomic ytterbium 
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We report on a measurement of a highly forbidden magnetic-dipole transition amplitude in 
ytterbium using the Stark-interference technique. This amplitude is important in interpret- 
ing a future parity nonconservation experiment that exploits the same transition. We find 
|(5d6s '^_Di|Ml|6s^ ^S'o)! = 1.33(6)stat(20)/3 x 10~*/.io, where the larger uncertainty comes from 
the previously measured vector transition polarizability /3. The Ml amplitude is small and should 
not limit the precision of the parity nonconservation experiment. 

PACS numbers: 32.70.Cs,32.60.+i,32.80.Ys 
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measure parity nonconservation 
-> 5c?6s ^Di transition in atomic yt- 



The proposal to 
(PNC) in theGs^iS'o 

terbium (Yb) |1| has prompted both theoretical H and 
experimental |4, ||] studies. The magnetic-dipole (Ml) 
amplitude for this transition is a key quantity for evaluat- 
ing the feasibility of a PNC-Stark interference experiment 
as proposed in Q. A nonzero Ml amplitude coupled with 
imperfections in the apparatus can lead to systematic un- 
certainties in a PNC experiment. Here we present the 
first experimental determination of the magnetic-dipole 
amplitude for the ^Sq — >^ Di transition. Our method is 
based on the technique of Stark interference H, 0, |§1 . 

In the absence of external fields, the ^5*0 Di tran- 
sition (Fig. |l|) is highly suppressed. An electric-dipole 
transition amplitude [A(i?l)] is forbidden except for par- 
ity mixing effects. A magnetic-dipole transition am- 
plitude is also forbidden because of the s-d nature of 
the transition. Consequently, a nonzero transition am- 
plitude exists only as a result of configuration mixing 
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FIG. 1: Low-lying energy levels of Yb. 



and spin-orbit interaction in both the upper and lower 
states Q. There have been no detailed calculations of 
this amplitude. Reference gives a rough estimate of 
|^(AfI)| < 10~*^oi where /Zq is the Bohr magneton. 

In the presence of a static external electric field, E, 
there is a parity-conserving mixing between the even 
parity ^Di state and the odd parity states. For a 
J — J — 1 transition, this mixing leads to a Stark- 
induced electric-dipole transition amplitude given by H 



A{E1 



St) 



i/3(E X f) 



(1) 



where e is the direction of the polarization of the laser 
light, (E X i)q is the q component of the vector in the 
spherical basis, and the vector transition polarizability /? 
is a real parameter. The magnitude of (3 was measured 



2.18(33) X 10-^eao/(V/cm). 



(2) 



In an electric field, the transition amplitude is the sum 
of the Stark-induced El amplitude and the forbidden All 
amplitude. The corresponding transition rate is 



W cx \A{Elst) + A{M1)\^ 

« \AiElst)\^ + 2Re[A{Elst)A{Ml)*] 



(3) 



where we neglect the contribution from |^(M1)|^ since 
|^(M1)| < |^(£;ist)| with the electric fields and polar- 
ization angles used here. The Stark-induced amplitude is 
proportional to E. Thus, reversing E changes the total 
transition rate, allowing the interference term to be iso- 
lated from the larger terms. The Ml amplitude is given 
by 



A{M1) = CDi, Mj\Ml\^So)(k X e)M.,, 



(4) 



where k is the direction of propagation of the excita- 
tion light. Equation |^ implies that only the Mj = ±1 
components of the upper state are excited by A{Elst), 
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where the axis of quantization is chosen along E. With 
k perpendicular to E, the sign of the interference term is 
opposite for the transitions to the Mj — ±1 components, 
as can be verified by a simple calculation. Thus, in or- 
der to observe the effect of the Stark-Ml interference, we 
apply a magnetic field, B, allowing us to resolve differ- 
ent magnetic sublevels. For B parallel to E (Fig. ^) the 
interference term in the transition probability is propor- 
tional to the rotational invariant (E x f) x (fc x e) • S. 

Comparison of the difference in the transition rate be- 
tween opposite electric field states to the sum is a mea- 
sure of the fractional asymmetry a, defined as 

W{E+) - W{E_) 

2(3 Ji, Mj|Ml|i^o)cos(g) ^,^ 
phi sin(f)) 

where is the angle between the dc electric field and the 
polarization of the excitation light (Fig. ||). The asym- 
metry increases with decreasing 6 while the dominant 
signal decreases as sin^ {9) . Most of the data was taken 
ai 9 = ±45°, where the interference term is maximal. 

Much of the apparatus used in this experiment had 
been used for the measurement of the Stark-induced tran- 
sition amplitude and is described in detail in Refs. |H, ||. 
A stainless steel oven with a multi-channel nozzle cre- 
ated an effusive beam of Yb atoms inside of a vacuum 
chamber with a residual pressure of w 5 x 10~^ Torr. 
The oven nozzle collimation resulted in a Doppler width 
for the 408-nm transition of ~ 150 MHz. The oven was 
heated with tantalum wire heaters operating at « 500 °C 
in the rear with the front « 100 °C hotter to avoid clog- 
ging. Ytterbium has seven stable isotopes with both zero 
and nonzero nuclear spin (^^^Yb, ^"^"Yb, ^''^Yb, ^'''^Yb, 
i^^Yb, / = 0; Yb, / = 1/2; and ^^^yb, / ^ 5/2). 
To avoid significant overlap of the optical spectra of the 
zero-nuclear-spin isotopes and the hyperfine components 
of the nonzero-nuclear-spin isotopes, an external vane 
collimator was installed; reducing the Doppler width to 
K, 15 MHz. The vane collimator was made by layer- 
ing 0.076 mm thick sheets of stainless steel foil between 
0.94 mm thick stainless steel spacers. The length of 
the collimator was 5.1 cm, providing a collimation an- 
gle of w 1°. The width of the collimator was 3.8 cm. 
The collimator was heated using tantalum wire heaters 
to ~ 350 °C to prevent clogging. The collimator was 
mounted on a movable platform, allowing precise align- 
ment of the angle of the collimator relative to the atomic 
beam during the experiment. We estimate an atomic 
density of « 2 x 10^ cm^"^ in the interaction region. 

Approximately 80 mW of laser light at 408 nm ex- 
cited ytterbium atoms to the ^Di state in the geometry 
shown in Fig. ^. The 408-nm light was produced by 
frequency doubling of w 1.2 W of 816-nm light from a 
titanium-sapphire laser (Coherent 899-21) pumped with 




FIG. 2: Schematic of apparatus. 



w 12 W from a multi-line argon-ion laser (Sectra Physics 
2080). A commercial bow-tie resonator with a Lithium- 
Triborate crystal (Laser Analytical Systems Wavetrain 
cw) provided frequency doubling. 

A Burle 8850 photomultiplier tube (PMT) monitored 
the fluorescence in the 6s6p ^Pi 6s^ ^Sq decay chan- 
nel at 556 nm (Fig. The top electrode had an array 
of 198 0.12 cm diameter holes, allowing the fluorescence 
to be collected by a Lucite light guide and conducted to 
the PMT. The presence of the holes in the top field plate 
reduced the electric field between the plates. This effect 
was calculated to be less than 1% using a random walk 
solution to Laplace's equation jl^. An interference filter 
with transmission centered at 560 nm with a 10 nm full 
width at half maximum was placed in front of the pho- 
tomultiplier tube in order to limit detection of scattered 
light at 408 nm. Approximately 0.05% of the atoms 
undergoing a transition were detected, resulting in typi- 
cal photocathode currents of lO'* e/s on the peak of the 
resonance for |E| = 15 kV/cm. The off-resonance back- 
ground was consistent with the expected PMT dark cur- 
rent and residual scattered light; contributing a noise w 3 
times less than the signal shot noise for |E| = 15 kV/cm. 
The laser frequency was scanned « 200 MHz with both 
increasing and decreasing frequency over the transition 
and the fluorescence spectrum was recorded with a dig- 
ital oscilloscope. The scan time each way was typically 
1 s. A typical single scan is shown in Fig. ^. 

After the laser was scanned the polarity of the electric 
field was either switched or left unchanged in accordance 
to the following pattern: (+ — — + — -I- -I- — ). 
This pattern was chosen to limit systematic effects asso- 
ciated with drifts in the laser frequency and atomic beam 
intensity. A bipolar power supply (Spellman CZEIOOOR), 
modified so that the polarity was computer controlled, 
produced the high voltage used in the experiment. The 
polarity of the top electrode was reversed while the bot- 
tom electrode remained grounded. A resistor divider 
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monitored electrode voltage. The magnitude of this volt- 
age changed by < 0.1% with the change in polarity. A 
delay of « 2 s after each switch allowed the electric field 
plates to fully charge before the next scan. The stain- 
less steel electric field plates were separated 1.016 cm by 
Delrin spacers. The typical value of the electric field was 
« 12 kV/cm. After each sequence of E-field switches, the 
polarity of the magnetic field was switched according to 
the (-1- — — h) pattern. The Earth's magnetic field was 
reduced to ^ 50 mG with external coils. A pair of in- 
vacuum coils in a near Helmholtz configuration provided 
the uniform magnetic field needed for the experiment. A 
typical magnetic field was 36 G. A run consisted of 5-10 
sets of 32 forward and backward laser scans (8 E-field 
switches x 4 B-field switches per set) with fixed values 
of 9, \E\, and \B\. Periodically « 50 laser scans were 
taken with zero magnetic field in order to monitor for 
changes in the lincshape of the transition due to temper- 
ature fluctuations of the oven and collimator. 

A temperature stabilized, hermetically sealed Fabry- 
Perot cavity with a free spectral range of 150 MHz and 
finesse of « 15 was used to monitor a portion of the 
816-nm light. A photodiode monitored the 408-nm laser 
power in order to normalize the signal for power fluctu- 
ations. The transmission through the Fabry-Perot and 
the 408-nm laser power were recorded concurrently with 
the fluorescence signal. 

The Fabry-Perot transmission peaks were used to line 
up the fluorescence spectra of two consecutive laser scans 
in order to eliminate frequency drift between scans. Two 
scans at opposite electric fields and the same magnetic 
field were combined. The sum of the two fluorescence 
spectra was fit to the function 

W{E+) + W{E^) = 

C(/(i^-Ci) + f{iy-C2))+miy + b, (6) 

where v is the frequency of the laser, Ci(2) is the center 
position of the first (second) peak, C is the amplitude of 
the peaks, m and b account for any linear background 
coming from scattered light. The function / was numeri- 
cally determined from the spectra taken at zero magnetic 
field. Because the sign of the interference term is oppo- 
site for the different magnetic sublevels (see Eq. ||), the 
spectral dependence of the asymmetry is given by the 
difference between the Mj — ±1 peaks multiplied by an 
asymmetry coefficient. The difference between the two 
spectra was therefore fit to a function whose line shape 
was constrained by the fit parameters from the sum fit: 

W{E+) ~ W{E-) = 

aaj{y-c^) - f{iy-C2)) + C, (7) 

where a is the asymmetry coefhcient given in Eq. ^, and 
C accounts, at lowest order, for any possible background 
that may be present in the difference due to a constant 
offset in the electric field which does not change sign with 
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FIG. 3: Typical power-normalized, single scan over the 
Mj = ± 1 components of the 6s^ ^5*0 ^ 5d6s "^Di transition 
of ^''■'Yb with |E| = 15 kV/cm, e = 45°, |B| = 36 G. 

the electric field switch. Higher order contributions to the 
line shape from a constant offset in the electric field were 
analyzed and found to be insignificant in the determina- 
tion of a. Changing the polarity of the magnetic field 
reverses the sign of a since the resonance frequencies of 
the magnetic sublevels switch (Eq. ||). Note that the 
asymmetry coefficient does not depend on the value of 
|B|. 

The measurements were performed on isotope ^^"^Yb 
which has a large relative abundance and is spectrally 
well isolated. The Ml transition amplitude was mea- 
sured in a variety of different field values and config- 
urations. The variation of E was from 5 kV/cm to 
20 kV/cm, e from -70° to 70°, and B from 12 G to 
84 G. In addition, data was taken without the exter- 
nal collimator in order to determine possible systematic 
effects associated with the line-shape modeling or fre- 
quency noise. For this data the overlap between the 
"4Yb and "3Yb (F = | ^ F' = |) fines was significant 
and the analysis modified to include effects of this tran- 
sition. 

The effects of misalignments of the fields and imper- 
fect reversals were analyzed analytically and numerically 
using density matrix formalism. These calculations indi- 
cate that systematic effects are significantly smaller than 
the statistical uncertainty. Possible systematic effects are 
also severely constrained by confirming the characteris- 
tics of the asymmetry. The method of analysis described 
above is sensitive to asymmetries which reverse sign with 
E and are of opposite sign for the two magnetic sublevels. 
Equation |^ implies that the sign of the asymmetry should 
also reverse with B and 9. Asymmetries which did not re- 
verse sign with either B or were consistent with zero. In 
addition, the dependence of the magnitude of the asym- 
metry on the magnitude of E and 9 was also verified. 

The final value of the Ml amplitude is based on data 
taken on two different days. The data is shown in Fig. |^. 
The statistical error for each point was estimated from 
the spread of values obtained for each complete sequence 
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FIG. 4: Experimental values of 1 2Mi | . The variation of E 
was from 5 kV/cm to 20 kV/cm and B from 12 G to 84 G. 
The circles (squares) represent data taken with (without) the 
external collimator. The solid line is the mean and the dashed 
lines are the statistical error on the mean. 



(« 100 times larger than in Cs |0|) makes the relative 
size of A{M1) to the PNC amplitude smaller than it is 
in Cs. An additional suppression of spurious interference 
between the Ml amplitude and the Stark-induced am- 
plitude is possible by using the geometry for the PNC 
experiment employed in Ref. ||l^. The reported mea- 
surement is for the isotopes with zero nuclear spin. The 
isotopes with nonzero nuclear spin (^^^Yb, 1=1/2 and 
lysy^j^ / = 5/2) have an additional contribution to the 
Ml amplitude and a small E2 amplitude due to hyper- 
fine mixing effects. However, these contributions are esti- 
mated to have values < IO^^/xq [0, and should only lead 
to small modifications to the present result. These effects 
will be investigated in future work. The size of the Ml 
amplitude should not limit the precision of a Yb PNC 
measurement which is in progress in our laboratory. 
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FIG. 5: Measured values of \^^\ as a function of |E|. The 
solid line is the expected dependence from the overall mean 
and the dashed lines correspond to the errors on the mean. 



of electric and magnetic field switches within a given con- 
figuration. These errors are consistent with the expected 
limit due to shot noise. Because the errors are estimated 
from the spread of 5-10 sets, there is some statistical vari- 
ation in the size of the error assigned to each run. There 
is additional variation of the errors due to differences in 
sensitivity for different polarization angles (see Eq. H) 
and differences in the amount of data taken in a given 
configuration. The final result is 



2{^D^, Mj = ±l\Ml\^So) 



-44.6(21)54,4 V/cm. (8) 



This corresponds to an Ml amplitude of 

|(3Di,Mj = ±l\Ml\'So)\ = 1.33(6)stat(20);3 x IQ-^ /lo, 

(9) 

where the second error represents the uncertainty in the 
determination of (3. 

The measured value of A{M1) agrees with the esti- 
mate in Ref. Q . This value is « 3 times larger than the 
corresponding amplitude in the cesium (Cs) transition 
where PNC is studied 
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|12| , |13| . However, the ex- 
pected large enhancement of the PNC amplitude in Yb 
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